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Abstract: We determine the optimal amount of life insurance for a household of two wage earners. 



We consider the simple case of exponential utility, thereby removing wealth as a factor in buying 
life insurance, while retaining the relationship among life insurance, income, and the probability of 
dying and thus losing that income. For insurance purchased via a single premium or premium payable 
continuously, we explicitly determine the optimal death benefit. We show that if the premium is 
determined to target a specific probability of loss per policy, then the rates of consumption are identical 
£S) • under single premium or continuously payable premium. Thus, not only is equivalence of consumption 

achieved for the households under the two premium schemes, it is also obtained for the insurance 
company in the sense of equivalence of loss probabilities. 
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1. Introduction 



The problem of finding the optimal strategy for purchasing life insurance has not been studied 
nearly as much as the corresponding problem for life annuities; see Milevsky and Young (2007) and 

in 

o\ . 

papers of Richard (1975) and Campbell (1980) and the more recent papers of Pliska and Ye (2007), 
l/^ ; Huang and Milevsky (2008), Kraft and Steffensen (2008), Nielsen and Steffensen (2008), Wang et al. 

(2010), Kwak, Shin, and Choi (2011), Bruhn and Steffensen (2011), and Egami and Iwaki (2011). We 
highlight the last two articles because that work is closely related to the work in this paper. Bruhn and 
Steffensen (2011) find the optimal insurance purchasing strategy to maximize the utility of consumption 
for a household, in which utility of consumption is a power function. Thus, the optimal death benefit is 
related to the wealth of the household. Similarly, Egami and Iwaki (2011) maximize household utility 
of consumption and terminal wealth at a fixed time and find that the optimal death benefit is a function 
of wealth because of the nature of their utility function. 

In this paper, we focus on the idea that, when maximizing the utility of household consumption, 
life insurance serves as income replacement; therefore, the amount purchased should be independent of 
the household's wealth. Indeed, when a wage earner dies, the wealth of the household remains, but the 
income of that wage earner is lost; life insurance can replace that lost income. Granted, life insurance 
is often awarded to the beneficiary in a single payment, but the beneficiary can use that payment to 
buy an annuity, which more clearly replaces lost income. 

To obtain an optimal amount of life insurance that is independent of the wealth of the household, 
we assume that utility of consumption is exponential. The resulting optimal death benefit is a function 
of the rates of mortality and rates of income of the household but not of its wealth, as desired. An 
added bonus of using exponential utility is that the optimization problem greatly simplifies, and we can 
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explicitly express the optimal death benefit, the optimal rates of consumption (both before and after 
the death of a wage earner), and the optimal investment in the risky asset. 

The rest of the paper is organized as follows: In Section 2, we model the household as two wage 
earners; the household wishes to maximize its (exponential) utility of consumption until the second 
death, with the possibility of purchasing life insurance that pays at the first death. Additionally, 
we allow the household to invest in a risky and a riskless asset. In Sections 3 and 4, we solve the 
optimization problem in the case of premium paid with a single premium or continuously until the 
first death, respectively. We find explicit expressions for the optimal amount of death benefit for the 
household to buy, and we determine how that death benefit varies with the underlying parameters. In 
Section 5, we compare the optimal rates of consumption from Sections 3 and 4 and show that they are 
identical if premium is determined by targeting a specific probability of loss per policy. 

For those readers who wish to skip the mathematical content of the paper, Section 6 provides a 
numerical example to illustrate our results, and Section 7 concludes the paper with a detailed summary 
of those results. We also have available an Excel spreadsheet at 

http://www.math.lsa.umich.edu/~erhan/OptimalLifeInsurance.xlsx 
that the interested reader can use to calculate the optimal amount of life insurance for his or her 
household. 

2. Financial Background, Statement of Problem, Verification Lemma 

In this section, we present the financial and insurance market for the household. Then, we state 
the optimization problem that this household faces. Finally, we present a verification lemma that we 
use to solve the optimization problem. 

2.1. Financial market 

The household invests in a Black-Scholes financial market with one riskless asset earning at the 
rate r > and one risky asset whose price process S = {S t }t>o follows geometric Brownian motion: 

dS t = fj,S t dt + a S t dB t , 

in which B = {Bt}t>o is a standard Brownian motion on a filtered probability space (0, J- , F = 
{•7"t}t>o> P)> with fx > r and a > 0. (As an aside, in courses on the theory of interest, actuaries call r 
the force of interest.) 

Let W t denote the wealth of the household at time t > 0. Let n t denote the dollar amount invested 
in the risky asset at time t > 0. Finally, let c t denote the continuous rate of consumption of the 
household at time t > 0. 

A consumption policy c = {c t }t>o is admissible if it is an F-progressively measurable process 
satisfying J* \c s \ ds < oo almost surely, for all t > 0. An investment policy II = {7r t } t > is admissible 
if it is an F-progressively measurable process satisfying f ir^ds < oo almost surely, for all t > 0. 
Furthermore, we require that the policies are such that Wt > — oo for all t > with probability 1. 
For example, this restriction prohibits the household from borrowing an infinite amount of wealth and 
consuming it all instantly. 

The household receives income from two sources: While family member (x) is alive, the household 
receives income at the continuous rate of I x > 0. Similarly, while family member (y) is alive, the 
household receives income at the continuous rate of I v > 0. Denote the future lifetime random variables 
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of (x) and (y) by t x and r y , respectively. We assume that t x and r y follow independent exponential 
distributions with means 1/X X and 1/X y , respectively. (In other words, [x) and (y) are subject to 
constant forces of mortality, X x and X y , respectively.) Let T\ denote the time of the first death of (x) 
and (y), that is, n = mm(T x ,T y ). Let r 2 denote the time of the second death, that is, r 2 = max(r x , r^). 

The household buys life insurance that pays at time n, the moment of the first death. This 
insurance acts as a replacement for the lost income. In this time-homogeneous model, we assume that 
a dollar death benefit payable at time t\ costs H at any time when both (x) and (y) are alive. Write 
the single premium as follows: 



H=(l + 6)A xy = (1 + 0) 



X x + X y 
X x + X y + r ' 



(2.1) 



in which 9 > is the proportional risk loading. Assume that is small enough so that H < 1 because 
if H > 1, then the buyer would be foolish to pay a dollar or more for each dollar of death benefit. The 
inequality H < 1 is equivalent to r > 9(X X + X y ). 

In this section and in Section 3, we suppose that the premium is payable at the moment of the 
contract; as stated above, H is the single premium. In Section 4, we consider the case for which the 
insurance premium is payable continuously at the rate h — (1 + — (1 + 0)(X X + X y ) until time t\. 

Let Dt denote the amount of death benefit payable at time t\ purchased at or before time t > 0. 
A life insurance purchasing strategy D = {D t } t > is admissible if it is an F-progressively measurable, 
non-decreasing process with W t > — oo for all f > with probability 1. Thus, with single-premium life 
insurance, wealth follows the dynamics 



dW t = (rW t - + {n~ r)7r t _ +I x +I y - ct_) dt + air t _ dB t — H dD t , 0<t<T 1 , 
W Tl =W Tl -+D Tl _, 

dW t = (rW t + (fj, - r)7r t + I x l{ Tl =r. a } + I y 1 {t 1 =t x } - c t ) dt + a-Kt dB t , n < t < r 2 . 



(2.2) 



A few words of explanation of the expressions in (2.2) are in order. Because we allow the life insurance 
purchasing strategy to jump due to lump-sum purchases, we write the subscript t— instead of t to denote 
the values of the corresponding process before any such jump. Before the first death, the household's 
wealth increases at the rate of interest earned on the riskless asset, r(W t _ — 7r t _), the drift on the 
risky asset, 7r t __ (modified by the random component cr7r t _ dB t ), and the rates of income, I x and I y . 
Also, before the first death, the household's wealth decreases with consumption and any purchase of 
life insurance. 

At the moment of the first death, the wealth of the household immediately increases by the death 
benefit D Tl _. After that death, the wealth process is much as it was before, except for the loss of the 
income of the family member who died. Also, the optimal amount of wealth invested in the risky asset 
and the rate of consumption chosen by the remaining family member will change. 

2.2. Utility of consumption 

We assume that the household seeks to maximize its expected discounted utility of consumption 
between now and r 2 , the time of the second death, by optimizing over admissible controls (c, II, D). 
The corresponding value function is given by 



U(w,D) 



sup E' 
(c,n,D) 



.D 



f 

Jo 



u(c t ) dt 



Ti > 



(2.3) 
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in which E 1 " 1 denotes conditional expectation given that Wo- = w and Do- = D. Here, r is the 
discount factor for the household, which we assume equals the rate of return for the riskless asset. The 
function u is the utility of consumption; we assume throughout this paper that 

u{c) = --e- ac , (2.4) 
a 

for some a > 0, which is called the coefficient of absolute risk aversion. Thus, exponential utility is 
equivalent to constant absolute risk aversion. 

To make this paper somewhat complete, we now present a way a household might determine its 
constant absolute risk aversion a. It is a method we borrow from Bowers et al. (1997). Suppose a 
household faces the possible loss of $10,000 with probability 0.01 and no loss with probability 0.99. In 
lieu of that random loss, the household is offered insurance for a fixed premium, and suppose P is the 
maximum dollar amount that the household would pay in exchange for the possibility of losing $10,000. 
That is, the household is indifferent between (i) paying P and receiving full reimburse for the loss and 
(ii) buying no insurance for the loss and paying for it out of pocket. Then, by expected utility theory, 
if the household's current wealth is w, the expected utility of these two alternatives are equal: 

e -a(w-P) = 0>01e -a(to-10000) + 0.QQ e -*w^ 

or equivalently 

P= - In (0.01e 10000a + 0.99) . 

For P £ (100, 10,000), this equation has a unique solution a > 0; see Gerber (1974) for more on the 
exponential premium principle. 

Remark 2.1. One can modify our set up to model a household composed of one income earner and 
other members. Indeed, let I x > and I y = 0. In that case, it would be reasonable to suppose that 
the household will buy insurance that pays only at time t x , instead of n . 

Additionally, one can modify our set up to model a family composed of two income earners and 
several children. In that case, let I x > and I y > 0. One could consider a variety of life insurance 
products in this case, for example, separate policies on (x) and (y). Also, it would make sense to 
maximize the utility of consumption until the third "death," namely, when the two wage earners and 
all the children have died. 

We leave the details of these adaptations to the interested reader and proceed with the case of a 
two-party household of two earners. We believe that the qualitative results we find here will hold in 
the cases suggested above. 

Remark 2.2. After the first death, the remaining family member, either (x) or (y), faces a Merton 
problem of maximizing utility of lifetime consumption. Specifically, the corresponding value function 
for the survivor, as a function of the hazard rate A and rate of income I of the survivor, is given by 

/■oo 

/ e-^^u^dt , 



in which wealth follows 



V(w;X,I) = sup W 

fc,ro 



/ 

Jo 



e rt u{c t )dt 



= sup E l 
(c,ro 



dW t = {rW t + r)ir t + I - c t ) dt + a n t dB t , W = w. 
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We know from Merton (1969, Section 9) that V is given by 

V(w; X,I) = — exp \ -ar ( w + - + — ^ J 1 , (2.5) 

ar { \ r ar 2 ) J 

in which m = \ {^-) ■ Note that wealth in equation (2.5) is augmented by the present value of a 
perpetuity that pays at the rate of J, in addition to a term that accounts for the force of mortality and 
the financial market. The optimal rate of consumption for the household after the first death is 

c* v (w;X,I) = ( u ')- 1 (V w (w)) = rw + I+^l, (2.6) 

a r 

which is the investment income if one were to invest all wealth in the riskless asset (rw), plus the 
survivor's earned income (J), plus an additional amount based on the survivor's force of mortality (A), 
the riskiness of the market (m), and the risk aversion (a). The optimal amount to invest in the risky 
asset after the first death is 

MW > A>/) = " ~ V^iwj = (2 - 7) 

a constant, which is the usual result when considering exponential utility. Not only is iXy a constant, 
it is independent of the force of mortality. 

Remark 2.3. It follows from the Dynamic Programming Principle (see, for example, Jeanblanc et al. 
(2004, Proposition 2.4)) that we can rewrite U in (2.3) as follows: 



U{w,D)= sup E 
(c,n,D) 



= sup E W ' D 
(c,n,D) 



l{r,<r y} QT" e- rt u(c t ) dt + e-^V(W Tx - + D Tx -;\y,I v )) 
+ 1 {r x > Ty} (J * e- rt u(c t )dt + e-^V(W Ty _+D Ty _;\ x ,I x )^J 

/ e- 5t [u{c t ) + \ x V{W t - + D t _;\y,I y ) + X y V{W t - + D t _; A x , /*)] dt 
Jq 



(2.8) 



in which 6 = r + A x + A y , u is given in (2.4), and V is given in (2.5). The first expectation includes the 
randomness in the times of death, while the second does not because the random times of death are 
accounted for within the integral. 

2.3. Verification lemma 

In this section, we provide a verification lemma that states that a smooth solution to a variational 
inequality associated with the maximization problem in (2.3) is the value function U. Therefore, we can 
reduce our problem to one of solving a variational inequality. We state the verification lemma without 
proof because its proof is similar to others in the literature; see, for example, Wang and Young (2012a, 
2012b). 

First, for c, n G R, define a differential operator C C,7T by its action on a test function / on R x R+; 
the definition of the operator is derived from the wealth dynamics before t\ in (2.2) and the alternative 
expression for U in (2.8). 



\<J 2 TT 2 f ww (w,D) - -, 
2 a 

(r + A x + Xy)f(w, D) + X x V(w + D; X y ,I y ) + X y V(w + D; X X ,I X 



f(w, D) = (rw + r)ir + I x + I y - c)f w (w, D) + ^a 2 ir 2 f ww (w, D) - -e^ 
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Lemma 2.1. Let T = T(w,D) be a function that is non- decreasing, concave, and twice- differentiable 
with respect to w and non- decreasing and differentiable with respect to D onRx R + . Suppose that T 
satisfies the following variational inequality: 

maxtmaxr^ T(w, D), T D (w, D) - H T w (w, D)) = 0. (2.10) 

C,7T 



Then, on R x R + , 

U = T. 

To specify the optimal strategies, first partition the region R x R + into TZ\ := {(w,D) G R x R + : 
U D (w,D)-HU w (w,D) < 0} anc2ft 2 := {(w,D) G R x R+ : U D (w, D) - H U w (w, D) = 0}. T/ien, i/ie 
optimal life insurance purchasing, consumption, and investment strategies are as follows: 
(a) When (W t _, -Dt_) G 7?-2, purchase additional life insurance AD so that (W t _ — HAD, D t _+AD) G 

(6) When (W t ,D t ) G cl(1Zi), purchase additional life insurance instantaneously to keep (W t ,D t ) G 
d(fti). 

(c) Consume continuously at the following rate when (W t , D t ) = {w, D) G cl(1Z\): 

c* u (w,D) = (u')- 1 (U w (w,D)). (2.11) 

(d) Invest the following amount of wealth in the risky asset when (W t ,D t ) = (w,D) G cl(Jl\): 

*( n\ l i ~ r U w (w,D) 

nuiw, D) = — - — — . (2.12) 

& U ww (w,D) 



The region TZ\ is called the continuation region because when the wealth and life insurance benefit 
lie in the interior of the region, the household does not purchase more life insurance; it continues with 
its current benefit. Indeed, Up < HU W means that the marginal benefit of buying more life insurance 
(Ud) is less than the marginal cost of doing so (HU W ). On the closure of that region, written cl{TZ\), 
the Hamilton- Jacobi-Bellman (HJB) equation holds: 

m&xL^ U(w,D) =0. (2.13) 

C,7T 

At the boundary of the continuation region, the household exercises singular control to remain 
in cl(1Z\). If the wealth and insurance benefit lie in 1Z 2 , then the household uses impulse control to 
move to the boundary of TZi, written d{Jl\). Because the household moves instantaneously from 1Z 2 to 
d(1Zi), we do not need to specify consumption and investment strategies on 1Z 2 . Also, because of this 
instantaneous movement, if (w,D) G 1Z 2 , then 

U(w, D) = U(w - H AD, D + AD), 

in which AD is such that (w — H AD, D + AD) G d(TZ\), so that the right-hand side is given by the 
solution to the HJB equation. 



3. Insurance Purchased with a Single Premium 
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Throughout this section, we assume that the household purchases its life insurance via single 
premiums H — (1 + 0) j-^j~+r P er dollar of insurance benefit, for some 6 > 0. In the next section, we 
consider the case for which the household pays its life insurance premium continuously until t\. 

3.1. Maximized utility 

On the closure of the continuation region, U = U(w, D) solves the HJB equation (2.13): 



(r + X x + X y )U = (rw + I x + I y )U w + max 



cU„ 



a 



+ max 



{fi - r)irU w + -a 2 ir 2 U 



,—ar(w+D)—m/r 



ar 



X x e -<*Iy-*v/r + \ e -al x -\ x /, 



(3-1) 

Suppose that U is of the form U(w,D) = e~ arw ; then, one can show k(D) is given implicitly as 

the unique positive solution of the following equation: 



k [rink + ar(I x + /„) + X x + \ y +m}= e -<*rD-m/r e -al v -\ v /r + Xy e - a I x -X x/r 

It will prove useful to know properties of k(D). 



(3.2) 



Lemma 3.1. k(D), defined as the unique positive solution of (3.2), is decreasing with respect to D, 
while k(D)e arD is increasing and convex with respect to D. 



Proof. To show that k'(D) decreases with respect to D, differentiate equation (3.2) to get 
k'(D) [r In k(D) + ar(I x + I v ) + X x + X y + m + r] 



-ar e 



-arD — m/r 



X x e - aI «- x y/ r + X v e - aI *~ x */ r 



(3.3) 



It follows from equation (3.2) that the expression in the square brackets on the left-hand side of (3.3) 
is positive; thus, k'(D) < 0. 

Next, from (3.2) and (3.3), observe that 



_d_ 

dD 



(k(D)e arD ) [r In k(D) + ar(I x + I y ) + X x + X y + m + r] = ar 2 k(D)e 



arD 



from which we deduce that the derivative of k(D)e arD with respect to D is positive. Differentiate this 
expression with respect to D to obtain 

d 2 d rk'(D) 

— (k(D)e arD ) [rlnk(D) + ar(I x + I y ) + X x + X y + m + r] + — (k(D)e^) ^L_2 

= ar 2 ±{k{D)e^ D ), 

which implies that the second derivative of k(D)e arD with respect to D is positive. □ 

Thus, the function (w,D) — > —^p-e~ arw is increasing, concave with respect to w and increasing 
with respect to D. According to Lemma 2.1, it is a good candidate for our maximized utility because 
it also solves the HJB equation (3.1). 

In the interior of the continuation region, we have Ud(w, D) — HU w (w, D) < 0, or equivalently, 



k(D) < exp 



H 
1-H 



a(I x + Iy 



X x + X y +m 



(3.4) 



8 



Note that the region defined by (3.4) is independent of wealth w, so in the case for exponential utility, 
the resulting control problem is not one of stochastic control. Either the household has enough life 
insurance D so that inequality (3.4) holds weakly, or the household buys enough insurance so that (3.4) 
holds with equality. 

In mathematical terms, if k(D) is greater than the expression on the right-hand side of inequality 
(3.4), then we are in the interior of the "buy" region IZ2, and the household will buy enough insurance 
(increase D) so that we have equality in (3.4). Recall that k'(D) < 0, so k(D) decreases as the household 
buys more insurance. 

We next find an expression for the boundary between the continuation and the buy regions, TZ\ and 
IZ2, respectively. In the case for exponential utility, this boundary is a line of the form {(w, D*) : w G R} 
for a fixed value of D* > 0. Despite the fact that the target death benefit, D* , is independent of wealth, 
it still depends on the parameters of the model in an interesting way. 

For D = D* , we have equality in (3.4), which gives us an explicit expression for k(D*). We, then, 
substitute this expression into (3.2) and solve for the optimal amount of life insurance D*: 



D* = max 



1 

ar 



In ^X x exp ^al x + + X y exp ^al y + ^ - In ( j^T^ 



1-H 







. (3.5) 



Note that we do not allow the household to buy a negative amount of life insurance; thus, we force 
D* > in the expression in (3.5). In the next proposition, we give the maximized utility in the case 
for which D* = in (3.5). In this case, the continuation region 1Z\ equals the entire domain R x R+. 

Proposition 3.2. If D* = in (3.5), then the maximized utility U defined by (2.3) on Rx R+ is given 
by: 

U(w,D) = -^-e- arw , (3.6) 
ar 

in which k(D) > uniquely solves (3.2). The associated optimal life insurance purchasing, consumption, 
and investment strategies are as follows: 

(a) Do not purchase more life insurance. 

(b) Consume continuously at the following rate when (W t ,D t ) = (w, D)eRx R + : 

c* ( W ,D) =rw- - lnk(D). (3.7) 
a 

(c) Invest the following amount of wealth in the risky asset when (W t ,D t ) = (w, D)eRx R + : 

' K u{ w i D) — ~~~2- (3-8) 
ar a z 



Proof. We use Lemma 2.1 to prove this proposition. Because U in (3.6) solves the HJB equation (3.1), 
to verify that equation (2.10) holds, it is enough to show that Ud{w, D) — H U w (w, D)<0onRx R+, 
which is equivalent to showing that (3.4) holding weakly. The left-hand side of (3.2) increases with 
respect to k; thus, to show that (3.4) holds, it is enough to show that the equality in (3.2) is replaced 
by > when we replace k with the right-hand side of (3.4). Specifically, it is enough to show that 

X x e -c<Iv-*v/r + \ e -aJ,-A,/rj ; 



rH 
1 - H 



exp 



H 



H 



a(I x + Iy) 



Ax + Xy 



g — arD — m/r 
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for all D G R + , or equivalently, 

exp (i~h) - K e " Ix+Xx/r + X v e a/y+A " /r , 



which holds if and only if D* = in (3.5). Thus, we have verified that U in equation (3.6) solves the 
variational inequality (2.10). Because the entire space is the continuation region, it is optimal not to 
purchase more life insurance. Finally, the optimal consumption and investment strategies are given in 
feedback form by the first-order necessary conditions of the HJB equation, as in equations (2.11) and 
(2.12), respectively. □ 

Next, we consider the case for which D* > in (3.5), and we state the main theorem of this section. 
Then, in Section 3.2, we examine how D* varies with the parameters of the model. 

Theorem 3.3. If D* > in (3.5), then the maximized utility U defined by (2.3) on Rx R + — TZ\WR,2, 
with Tlx = {(w, D) : w G R, D > D*} and TZ 2 = {(it), D) : w G R, < D < D*}, is given by: 

(i) For (it), D) G d{Ki) = {(w, D) : w G R, D> D*}, 



U(w,D) = -^-e- arw , (3.9) 
ar 



in which k(D) > uniquely solves (3.2). 
(it) For (w,D) G TZ 2 , 

U(w, D) = U(w - H(D* - D), D*). (3.10) 

The associated optimal life insurance purchasing, consumption, and investment strategies are 
as follows: 

(a) Purchase immediately additional life insurance of AD = D* — D when D t - = D < D* . 

(b) Do not purchase additional life insurance when D t - = D > D* . 

(c) Consume continuously at the rate given by (3.7) when (Wt,Dt) = (w,D) G ci(TZ\). 

(d) Invest the amount of wealth in the risky asset given by (3.8) when (Wt, D t ) = (it), D) G cl(TZ\). 

Proof. As in the proof of Proposition 3.2, we use Lemma 2.1 to prove this theorem. By construction, 
U in (3.9) satisfies the HJB equation on cl(H\). We next show that U in (3.10) satisfies C C,1T U < 
on TZ 2 for all c, n G R. To this end, observe that for (w,D) G TZ 2 , U w (w,D) = U w (w - HAD,D*) 
and similarly for the second derivative with respect to w. Thus, for (w,D) G TZ 2 , we have from the 
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definition of C c ' n in (2.9), 

£ c ' n U(w, D) = (rw + (ti- r)n + I x + I y - c)U w (w, D) + ^a\ 2 U ww (w, D) + u(c) 

- (r + X x + X y )U(w, D) + X x V(w + D; X y , I y ) + X y V(w + D; X x , I x ) 
= (r(w - H AD) + (fi- r)ir + I x + I y - c)U w (w, D) + ^a 2 n 2 U ww (w, D) + u(c) 

- (r + X x + X y )U(w, D) + X x V(w + D*;X y ,I y ) + X y V(w + D*;X X ,I X ) 

+ rH AD U w (w, D) + X x (V(w + D;X y , I y ) - V(w + D*;X y ,I y )) 

+ X y (V(w + D; X X ,I X ) - V(w + D*;X X ,I X )) 
= C c ' n U(w - H AD, D*) + rH AD U w (w, D) 

+ X x (V(w + D; X y ,I y ) - V(w + D*;X y , I y )) + X y (V(w + D; X X ,I X ) - V(w + D*; X X ,I X )) 
< rH AD U w (w-H AD, D*) + X x (V(w + D; X y ,I y ) - V(w + D*; X y , I y )) 

+ X y (V(w + D; X X ,I X ) - V(w + D*;X X ,I X )) 
= rH AD exp { ^ - a(I x + I y ) - X * + X v + m | e -«r(«,-fl ad) 

1 / TO S 

ar 



exp (-arw - ^) (e- arD - e - arD ^ X x exp (-al y - + X y exp (-al x - 



in which the inequality follows from the fact that £ c,7r U < on cl(1Z\). Also, in the last expression, 
we substitute for both U and V and use the fact that k(D*) for D* > is given by the right-hand side 
of (3.4). After simplifying, we learn that the last expression is non-positive if and only if 

rH AD e ^HAD +T M w < J_ ^- arD _ e - arD > ^ K exp ^ + ^ +XyGxp ^ aIy + ^ 



From the definition of D* > in (3.5), substitute for X x exp (al x + ^f) + X y exp (al y + ^ E 
above to get the following inequality: 



in the 



riJADe arMD+ A < — ( e - arD 

ar V 

or equivalently, 



-arZT\ ocrD* rH 

ar V J 1 - H 



e -arAD + ^ _ A£) e -or(l-fl)AD < ^ 

in which we rely on the assumption that H < 1. This last inequality holds if e~ a + abe~ ab < 1 for any 
a > and b £ (0, 1), which, via elementary analysis, one can show is true. Thus, we have demonstrated 
that max Cj7r £ c,7r U < on R x R + with equality on cl(lZi). 

Next, we show that Up — H U w < on R x R + . By construction, this inequality holds with 
equality on 1Z 2 ; see equation (3.10). Now, XJ D — H U w < if and only if inequality (3.4) holds. By 
definition of D* , k(D*) is given by the right-hand side of (3.4). Recall that k(D) strictly decreases with 
respect to D; thus, inequality (3.4) holds exactly when D > D* , or equivalently on TZi. Thus, we have 
shown that Ud ~ HU W < on R x R + with strict inequality on 1Z\ and equality on TZ 2 - 

We have verified that U given in equations (3.9) and (3.10) satisfies the variational inequality 
(2.10) in Lemma 2.1. Also, U is non-decreasing, concave, and twice-differentiable with respect to w and 
non-decreasing and differentiable with respect to D. Therefore, U as stated is the maximized utility, 
and the optimal strategies follow from Lemma 2.1. □ 
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Remark 3.1. Because the optimal death benefit D* is independent of wealth, one can alternatively 
determine D* as follows: First, suppose the household begins with wealth w and no life insurance. If 
the household then buys life insurance of amount D, the maximum utility equals U(w — HD,D) in 
which U is given by equation (3.6) or (3.9). Next, find D = D* > to maximize U(w — HD,D), or 
equivalently, find D* > to minimize k(D)e arHD . Because k(D)e arHD is convex with respect to D, 
there is a unique value of D* > that minimizes it, and one can show that that value is the one given 
in (3.5). 

We have the following corollary to Proposition 3.2 and Theorem 3.3, in which we compare the 
optimal consumption and investment immediately before and after the first death. 

Corollary 3.4. At the moment of the first death, the optimal rate of consumption changes by 

Ac(D)=rD + I a + ^^ + -lnk(D), (3.11) 

ar a 

in which I a and X a are the income rate and force of mortality of the survivor, respectively. When 
D = D* > 0, then the expression in (3.11) becomes 



Ac(D*) = - 

a 



A d Aa ( ,j M , A„-A d \\ fr-6(\ x + \ v ) 
In - — — — + - — — — exp a(I a - I d ) H + In 1 



Az + Ay \ x + \ y V r 77 V r(l + 9) 



(3.12) 



in which the subscripts a and d refer to the one who is alive or dead, respectively. The optimal amount 
invested in the risky asset does not change upon the first death. 

Proof. From Proposition 3.2 and Theorem 3.3, we know that the optimal rate of consumption immedi- 
ately before the first death equals rw — ^ In k(D). From equation (2.6), the optimal rate of consumption 
immediately after the first death equals r(w + D) + I + (A + m)/(ar), in which I and A are the in- 
come rate and force of mortality of the survivor, respectively. The change in the rate of consumption, 
therefore, equals the expression in (3.11). 

When D* > 0, the right-hand side of inequality (3.4) equals k(D*) and (3.5) gives us D*. After 
substituting these expressions into (3.11), we obtain (3.12). From equations (2.7) and (3.8), we see that 
the optimal amount invested in the risky asset is constant, before and after the first death, independent 
of wealth and the death benefit. □ 

Numerical work indicates that the change in consumption can be positive or negative. We deduce 
the following properties of the change in consumption from the expression in (3.11): 

(i) If the survivor has greater income and greater force of mortality than the household member who 
died, then the change in consumption is greater than if the latter had survived. This result makes 
sense because the surviving household member has more income from which to consume and has 
a shorter expected time until dying so will also consume more for that reason. 

(ii) From (3.2) and (3.3), one can show that arD + lnk(D) increases with D; thus, Ac(D) increases 
with D. So, for a household with a larger death benefit, that household's change in consumption 
will be greater than for a household with a smaller death benefit. This result, too, makes sense 
because-all else equal-one expects that a household receiving a larger death benefit will consume 
more. 

(iii) As a special case of observation (ii), if a household does not buy any life insurance, then the change 
in the rate of consumption at the first death is Ac(0), which is less than Ac(D) for D > 0. 
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3.2. Properties of the optimal death benefit 

In this section, we determine how D* varies with the parameters of the model. Most results are as 
one might expect; however, one of them is rather surprising. Indeed, in Corollary 3.7, we find that if the 
household is close to risk neutral but still risk averse (a > but small), then buying no life insurance 
is optimal, even when it is priced fairly (6 = 0). We also learn in Proposition 3.10 that the optimal 
amount of insurance is bounded above by the maximal loss, as measured by the present value of the 
perpetuity that pays at the rate of m&x(I x , I y ), an intuitively pleasing result. 

First, note that although the optimal consumption and investment strategies depend on the drift 
and volatility of the risky asset, the optimal death benefit D* does not. 

Proposition 3.5. The optimal amount of life insurance D* decreases with the premium loading 6 to 
the extent that for 6 large enough, buying no life insurance is optimal. 

Proof. From the expression for D* in (3.5), define / by 

rH{9) \ H{9) 



/(<?) = - In 



i-H(6)j \-H{ey 

in which H = H(9) is given in (2.1). It is enough to show that / decreases with 6. 

J{> H(6){l-H{6)Y ' 

because H'(9) > 0. From (3.5), we see that as H approaches 1, D* goes to 0. □ 

It is no surprise that the optimal amount death benefit decreases as the loading on the single 
premium increases. Proposition 3.5 helps to confirm the reasonableness of our model's qualitative 
results. Similarly, we expect that the optimal death benefit increases as the household becomes more 
risk averse, and this is indeed the case, which we prove in the next proposition. 

Proposition 3.6. The optimal amount of life insurance D* increases with the risk aversion of the 
household, as measured by a. 

Proof. From (3.5), define g by 

1 



a 



In ( A, exp ( al x + j + X y exp (al y + J J - In (-^Jj ) ~ y^jj 



(3.13) 



then, 



g'(a) = -i 



+ 



In ^A x exp (al x + ^ + X y exp (al y + - In (f^) " Y~H 

1 I X X X exp (al x + + I y \y exp (al y + 



a X x exp (al x + ^f) + X y exp (al y + ^ 
which is positive if and only if 

In [X x exp (al x + — ) + X y exp ( al y + — ) ) - In ' 



H 1-H 



< a 



I X X X exp (al x + ^) + I y X y exp (al y + ( 3 - 14 ) 
X x exp (al x + ^)+X y exp (al y + ^ 
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When a = 0, one can show that the left-hand side of (3.14) is negative. Indeed, when 8 = and a = 0, 
it reduces to In ^ A ^ A e~ A "/ r + A X + x e~ A */ r j < 0. Because the left-hand side of (3.14) decreases with 
respect to 8 when H < 1, it is negative for all 8 > such that H < 1. 

Because inequality (3.14) holds when a = 0, it is enough to show that the derivative of the left- 
hand side with respect to a is less than or equal to the derivative of the right-hand side. This ordering 
of the derivatives is equivalent to 



>0, 



d I X X X exp [al x + If) + IyX y exp ( al y 



da 



X x exp (al x + ^)+Xy exp (al y + ^ 



which is straightforward to demonstrate. □ 

From the proof of Proposition 3.6, we have the following corollary. 

Corollary 3.7. For 8 > 0, there exists a(8) > such that purchasing no insurance is optimal if and 
only if a < a(8). 

In other words, even for actuarially fair insurance (that is, 8 — 0), if the household is not very 
risk averse, then it will not buy life insurance. This result contrasts dramatically with the one in, say, 
casualty insurance in which a risk averse agent will buy full coverage for a loss if insurance is priced 
fairly; see, for example, Exercise 1.22 in Bowers et al. (1997). 

In the next proposition, we explore when the optimal amount of life insurance is concave with 
respect to a. That is, as a increases, the optimal death benefit increases (as we know from Proposition 
3.6) but at an increasingly slower rate. 

Proposition 3.8. If (I X ,X X ) > (I y ,X y ) or if (I X ,X X ) < (I y ,X y ), then the optimal amount of life 
insurance D* is concave with respect to the risk aversion of the household, as measured by a, when 
D* > 0. 

Proof. It is enough to show that g in (3.13) has a negative second derivative with respect to a. 



5» = A 



, ( T KW , ( rH \ H 
In X x exp al x H + exp al y + — - In 



r ) y " V r ) ) \1-H J l-H 

2 I X X X exp (al x + ^) + ly^y ex P ( a ^y + T") 
* 2 X x exp(al x + ^) +X y exp(aI y + ^) 

(4 - I y ) 2 X x X y exp (a(I x + I y ) + 



1 



' X x exp (al x + ^) + X y exp (al 

which is negative if and only if 

, ( , ( r A x \ x f Xy\\ ( rH \ H 
In A^ exp al x H + X y exp al y + — - In 



< a 



r J y "V r ) ) \l-H J l-H 

I X X X exp (al x + ^) + I y Xy exp (al y + ^) a 2 (I x - I y ) 2 X x X y exp (a(I x + I y ) + ^±^) 



X x exp(aI x + ^)+X y exp(aI y + ^ 2 (\ x exp ( a / x + >f) + X y exp (al y + ^ 



(3.15) 
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As we observed in the proof of Proposition 3.6, the left-hand side is negative when a = 0; thus, it 
is enough to show that the derivatives with respect to a of the two sides of (3.15) are ordered. This 
ordering of the derivatives is equivalent to 



d (Ix-Iy) 2 K\exp(®(I x + Iy) + 

da [x x exp (al x + ^)+X y exp (al y + ^ 



< 0, 



which reduces to 



(Ix - ly) ( ^ ex P ( aI x + — 



X y exp ^al y + — ^ ^ > 0. 



This last inequality holds if either (I x , X x ) > (I y , X y ) or if (I x , X x ) < (I y , X y ), and we are done. □ 

Even though ordering of the derivatives in (3.15) is a stronger condition than one needs to prove 
that g is concave with respect to a, numerical experiments indicate that g might not be concave when 
the conditions of Proposition 3.8 do not hold. 

Next, we observe that D* is increasing and convex with respect to the incomes of the household; 
we state this result without proof because it is clear from the expression for D* in (3.5). 

Proposition 3.9. When D* > 0, the optimal amount of life insurance is increasing and convex with 
respect to I x and I y . 

The sign of the derivative of D* with respect to either force of mortality is ambiguous. Indeed, 
see the expression for D* in (3.5): increasing X x or X y shortens the expected horizon, which tends to 
increase D* via the logarithmic term; on the other hand, increasing either force of mortality increases 
the cost of life insurance, which tends to decrease D* via the i7-terms. 

We end this section with an upper bound on the life insurance purchased by a household. 

Proposition 3.10. 

max(/ x , Jj,) 
D < u — = max(4, Iy) aoo| • 

That is, D* is bounded above by the present value of a perpetuity that pays at the rate of max(/ x , I y ). 
Moreover, linic^oo D* = max(/ x , I y ) . 

Proof. From Proposition 3.5, we know that D* for 9 > is bounded above by the D* corresponding 
to 9 = 0. Thus, from (3.5), we have 



D* < max 



< max 



1 

ar 
1 

ar 



In 



A, 



Ax + Xy 

a m&x(I x ,I y ) 



exp ( al x + — ) + 



A, 



A X ' + Ay 



exp I al y H — - 



A x - + A, 







In 



A., 



Ax + A^ 



-X y /r 



+ 



A, 



Ax + A^ 



1 / / -r -r n niax(/x> Iy) 

< — (a max = V y > . 

ar r 

By applying L'Hopital's rule to the expression for D* in (3.5), we obtain lmic^oo D* = m&x(I x , I y ) ex- 



it is interesting that the optimal amount of death benefit is bounded above by the maximal "loss" 
that the household might sustain, in which we view the loss as the present value of the income stream. 
Thus, it is never optimal for the household to buy "too much" insurance. 
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4. Insurance Purchased by a Continuously Paid Premium 

In this section, we modify the work of Sections 2 and 3 by assuming that the household buys its 
life insurance via a premium paid continuously until the first death at the rate of h — (1 + 0)(X X + X y ) 
per dollar of insurance for some 6 > 0. 



4-1. Maximized utility 

As in Section 3.1, the optimal strategy for purchasing life insurance is for the household to buy 
additional insurance of D - D if D < D for some optimal amount of insurance D ; otherwise, 
the household buys no further insurance. In this section, we demonstrate this result and provide an 
expression for D* . 

With continuously paid premium for life insurance, wealth follows the dynamics 

(- dW t = [rW t - + (ji - r)7r t _ + 4 + I y - ct- - hD t -) dt + a 7r t _ dB t , < t < n, 

W T1 =W T1 _+D T1 _, (4.1) 
, dW t = (rW t + (n - r)7r t + I x l{ Tl=T , a } + I y l{ n =r x } ~ <h) dt + an t dB t , n<t< r 2 . 

By comparison with the wealth process in (2.2), note that the only change is to replace the single- 
premium term —HdDt with the continuously-paid-premium term —hDt-dt. 

Denote the maximized utility function by U, which is defined as in equation (2.3) with the wealth 
dynamics given in (4.1). Its HJB equation when D > D* is similar to the one for U given in (3.1) with 
an additional term on the right-hand side equal to —hDU w . Thus, on R x R+, U solves the variational 
inequality: 

max(max£ c ' 7r U(w,D) - hDU w (w,D), U D {w,D)) = 0. (4.2) 

C,7T 

The analog of the verification lemma, Lemma 2.1, holds in this case. 

If we hypothesize a solution of the form U(w,D) = — ^p-e~ arw , then k(D) is given implicitly as 
the unique positive solution of the following equation: 

k [rink - arhD + ar(I x + Iy) + X x + X y + m] = e -^D- m /r N e - a i v -\jr + Xy e ~ai x -\ x /A ^ 

Note equation (4.3) 's similarity to (3.2). 

On the continuation region, we have U d < because on that region it is optimal not to purchase 
life insurance; that is, by increasing D, we would decrease U. Note that Ud < is equivalent to 
k'(D) > 0, which one can show occurs if and only if D > D* , in which D* is given by 



D — max 



a(h + ■ 



In ^X x exp ^al x + + X y exp ^al y + ^ - In 



K- h - 

r 







(4.4) 



the analog of D* in (3.5). In the case for which 6 = = 6, we have j-^j = 7, so in that case, the only 
difference between (3.5) and (4.4) is the divisor of ar versus a(h + r) = a(X x + X y + r), respectively. 
Note that 1/r = a^, while 1/(X X + X y + r) = d xy . 

As in Section 3.1, we have two cases to consider: (1) D* = 0, and (2) D* > 0. The next two results 
are the analogs of Proposition 3.2 and Theorem 3.3, respectively. Their proofs are similar to the ones 
in Section 3.1, so we omit them. 
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Proposition 4.1. If D = in (4.4), then the maximized utility of lifetime consumption U on Rx R+ 
is given by: 

U(w,D) = -^-e- arw , (4.5) 
ar 

in which k(D) > uniquely solves (4.3). The associated optimal life insurance purchasing, consumption, 
and investment strategies are as follows: 

(a) Do not purchase more life insurance. 

(b) Consume continuously at the following rate when (W t ,D t ) = (w, D)eRx R + : 



dUw,D) = rw--]nk(D). (4.6) 

u a 



(c) Invest the following amount of wealth in the risky asset when (W t ,D t ) = (w, D)eRx R + : 

^w,D) = t^L. (4.7) 
u ar a 2 



Theorem 4.2. If D > in (4.4), then the maximized utility of consumption U on R x R+ = 1Z\ UIZ2, 
with Tlx = {{w,D) :w£~R,D >D*} andTZ 2 = {{w,D) : w G R, < D <D*}, is given by: 
(i) For (it), D) G cl(TZi) = {(w, D) : w G R, D> D* }, 

U(w,D) = -^l e - arw , (4.8) 
ar 

in which k(D) > uniquely solves (4.3). 
(it) For (w,D) G K 2 , 

U(w,D) = U(w,D*). (4.9) 

The associated optimal life insurance purchasing, consumption, and investment strategies are 
as follows: 

(a) Purchase additional life insurance of AD = D* — D when D t _ = D < D* . 

(b) Do not purchase additional life insurance when D t _ = D > D* . 

(c) Consume continuously at the rate given by (4.6) when (W t ,D t ) = (w,D) G cl(TZ\). 

(d) Invest the amount of wealth in the risky asset given by (4.7) when (Wt, D t ) = (w, D) G cl{TZ\). 

Remark 4.1. As in the case of a single premium, there is an alternative way to find the optimal death 
benefit D* ; see Remark 3.1. Because the optimal death benefit is independent of wealth, one can find 
D as follows: First, for a given death benefit D, the maximum utility equals U(w, D) in which U is 
given by equation (4.5) or (4.8). Then, find D = D > to maximize U(w,D), or equivalently, find 
D > to minimize k(D). Because k(D) is convex, there is a unique value of D > that minimizes 
it, and one can show that that value is the one given in (4.4). 

Parallel to Corollary 3.4, we have the following corollary to Proposition 4.2 and Theorem 4.2, 
in which we compare the optimal consumption and investment immediately before and after the first 
death. We also compare this change with the corresponding change when premium is paid in a lump 
sum for a given death benefit. 
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Corollary 4.3. At the moment of the first death, the optimal rate of consumption changes by 

Ac(D) = r D + I a + + - lnk(D), 

ar a 

in which I a and X a are the income rate and force of mortality of the survivor, respectively. As before, 
the optimal amount invested in the risky asset does not change upon the first death. Furthermore, for 
a given death benefit D > 0, Ac(D) > Ac(D) with equality only when D = 0. 

Proof. We rewrite equations (3.2) and (4.3), respectively: 

rk{D)e a ^+ 1 ^ b£± ^ 1 In (k{D)e^+ I ^ h ^^) = e~ arD (a.,,"' ~ + \ y e a1 ^) , 

and 

rk(D)e<*( I * + ^ +h * ± ^ [in (k(D)e a ^ +I ^ h ^ ± ^) - ahD] = e^ D (a.,,"' ^ + \ y e^+^) . 

Thus, k(D) > k(D) with equality only when D = 0, which completes our proof. □ 

That the rate of change of consumption is greater in the case of continuously paid premium is 
not a surprise because part of the outflow of this household before the time of the first death is the 
premium, so that the first death not only triggers receipt of the death benefit D, it triggers the end of 
the premium payment. 

4-2. Properties of the optimal death benefit 

In this section, we determine how D varies with the parameters of the model. Most of the 
results of Section 3.2 hold in this case because the expression for D is similar to one for D* , including 
the surprising Corollary 3.7, in which we found that a risk-averse household might not buy any life 
insurance, even when it is priced actuarially fairly. We state the following propositions and corollary 
without proof. 

As in the case of single-premium life insurance, although the optimal consumption and investment 
strategies depend on the drift and volatility of the risky asset, the optimal death benefit D does not. 

Proposition 4.4. The optimal amount of life insurance D decreases with the premium loading 9 to 
the extent that for 9 large enough, buying no life insurance is optimal. 

It is no surprise that the optimal amount death benefit decreases as the loading on the continuous 
premium increases. Similarly, we expect that the optimal death benefit increases as the household 
becomes more risk averse, and this is indeed the case, which we state in the next proposition. 

Proposition 4.5. The optimal amount of life insurance D* increases with the risk aversion of the 
household, as measured by a. 

As in Section 3.2, we have the following corollary to Proposition 4.5. 

Corollary 4.6. For 9 > 0, there exists a(9) > such that purchasing no insurance is optimal if and 
only if a < a(8). 

In other words, even for actuarially fair insurance (that is, 9 = 0), if the household is not very risk 
averse, then it will not buy life insurance, a result that we did not expect in this case nor in the case 
of single-premium life insurance. 
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In the next proposition, we find that the same conditions as stated in Proposition 3.8 ensure that 
the optimal amount of life insurance is concave with respect to a. 

Proposition 4.7. If (I X ,X X ) > (I y ,X y ) or if (I X ,X X ) < (I y ,X y ), then the optimal amount of life 
insurance D* is concave with respect to the risk aversion of the household, as measured by a, when 
D* > 0. 

Next, we observe that D > is increasing and convex with respect to the incomes of the household; 
we state this result without proof because it is clear from the expression for D in (4.4). 

Proposition 4.8. When D > 0, the optimal amount of life insurance is increasing and convex with 
respect to I x and I y . 

The sign of the derivative of D* with respect to either force of mortality is ambiguous. Indeed, 
see the expression for D* in (3.5): increasing X x or X y shortens the expected horizon, which tends to 
increase D* via the logarithmic term; on the other hand, increasing either force of mortality increases 
the cost of life insurance, which tends to decrease D* via the i7-terms. 

We end this section with an upper bound on the life insurance purchased by a household. 

Proposition 4.9. 

D < ^f- r y > < max(/ x , I y ) a xy . 

That is, D* is bounded above by the present value of a joint life annuity on (xy) that pays at the rate 
o/mas(4,/ y ). Moreover, lim^oo D* = max(4,/ 9 ) a xy . 

Thus, we see that an upper bound for D* is the actuarial present value of a joint life annuity 
that pays at the rate of m&x(I x , I y ); this contrasts with the result of Proposition 3.10, namely, D* < 
m&x(I x ,I y )a^. 

5. Actuarial Considerations 

In this section, we consider two problems inspired by actuarial mathematics. The first is to compare 
the rates of consumption when the single premium and continuous premium are calculated to ensure 
a given probability of loss per policy. The second problem is to find the probability that consumption 
reaches zero, which relates to the hitting time of Brownian motion with drift, a basic problem in ruin 
theory. 

5.1. Targeting a Specific Probability of Loss 

For a given household, we compare the rates of consumption when the single premium H and 
continuously-paid premium h are determined by a given probability of loss. Assume that the household 
begins with Wo- = w wealth and Do- = death benefit and that the optimal amounts of death benefit 
for both cases are positive, that is, D* > and D* > 0. 

In the case of a single premium, the household immediately spends HD* so that Wo = w — HD* 
and Dq = D*. In fact, D t = D* for all t > 0. By following the wealth dynamics in (2.2) for the 
optimally controlled wealth, as stated in Theorem 3.3 before the first death and as given in equations 
(2.6) and (2.7) after the first death, we obtain the optimally controlled wealth for premium payable as 
a lump sum: 

( w -HD* + 6t + ^B t , t<n. 
1 \w* i + ^(t-n) + ^(B t -B Tl ), n <t<n, ^ ) 
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in which 

2m 1 
5 = — +I x +I y + -hik(D*). (5.2) 

ar a 

The corresponding optimal rate of consumption is given by 

f r(w - HD*) + rSt + ^B t - ± In k(D*), t<n, 
1 \r{w-HD*)+rD* +r5r 1 + ^{t-T l ) + ^B t +I a + ^, n < t < r 2 . 

Similarly, the optimally controlled wealth for premium payable continuously until the first death 
equals 

( w + St + ^B t , t< n . 

* In+^i-Til + ^^-fl,), r 1 <t<75, l ' J 

in which 

2777 I J; * 

<5 = — + 4 + I v + - In k(D )-hD . (5.5) 
ar a 

The corresponding optimal rate of consumption is given by 



rw + r5t + ^B t -±lnk(D ), t<n, 

rw + rD* + rdn + - n) + ^B t + I a + **±E», n < t < r 2 . 



(5.6) 



Next, we suppose that the single premium H and the premium h are chosen by the insurance 
company to target a specific probability of loss q on each policy. The loss-at-issue random variables are 
L = D*e~ rTl — HD* and L = D*e~ rTl — hD*a—\, respectively. Because T\ is exponentially distributed 
with mean \ j(\ x + A y ), it is straightforward to show that q = P(L > 0) implies that 



q= l-H—^- 



or equivalently, 

H = (l-q) T ^j. (5.7) 
To ensure that H > ' we rec l u i re that Q < 1 ~~ (; 



_ \ x +X v +r 

Similarly, one can show that q = P(L > 0) implies that 



q = l 

or equivalently, 



h 



h + r 



r(l-q)^+^y 



l-(l-g) 



(A -hA \ r 
a +a +r ) ' wn i cn i s the identical condition 

for > . A lt • Therefore, throughout the rest of this discussion, we assume that 



'^-(XT^I ■ (5 - 9) 



The proof of the next lemma is easy, so we omit it. 
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Lemma 5.1. // H and h are given by equations (5.7) and (5.8), respectively, then h = j^L an d 
H = -rh- 

h+r 

The next lemma we be useful in proving our main result concerning c* t and cj* . 
Lemma 5.2. If H and h are given by equations (5.7) and (5.8), respectively, then 5 = 5. 
Proof. From the right-hand side of inequality (3.4), we have 

lnfc ^*) = J~fj ~ a ^ + *v) ~ Xx + X r y + m , (5.10) 

and similarly, 

\nk(D*) = - r + ahD* - a(I x + I y ) - A * + + m (5 . n) 
Thus, by substituting (5.10) and (5.11) into (5.2) and (5.5), respectively, we obtain 

S= ™-**-*v + LH (5.12) 
ar a 1 - H 

and 

_^m-\ x -X y + lh 
ar ar 

It follows from Lemma 5.1 that 5 = 5. □ 

We give the main result of this section in the following theorem. 

Theorem 5.3. If H and h are given by equations (5.7) and (5.8), respectively, then c* = ~c[ for all 
<t < t 2 with probability 1. 

Proof. From equations (5.3) and (5.6), we learn that c* = cj. for all t < t\ if and only if rHD* + 
£ lnk(D*) = ± In k(D ), in which we use the fact that 5 = 5. This condition reduces to rHD* = hD , 
which follows from equations (3.5) and (4.4) by using h = j^- and H = j^-p from Lemma 5.1. Thus, 
c* t = c* t for all t < ti . 

We proceed similarly to show c* = c\ for all T\ < t < r 2 ; this equality is equivalent to (1 — H)D* = 
D* , which follows from equations (3.5) and (4.4) by using Lemma 5.1. We have, thus, proved the 
theorem. □ 

From the proof of the above theorem, we note that D* > D* ; however, because the rates of 
consumption are identical under either premium-payment scheme, the household would be indifferent 
between the two. We have the following corollary to Theorem 5.3. 

Corollary 5.4. If 9 = 8 = 0, then c* = ~c% for all < t < r 2 with probability 1. 

Proof. The condition that 8 = 8 = is equivalent to H and h determined by targeting the probability 
of loss equal to the right-hand side of inequality (5.9). □ 

Remark 5.1. We find it quite surprising that the two rates of consumption are identical when the single 
premium and premium rate are computed by a loss-percentile principle. One might be tempted to say 
that two risk loadings 8 and 8 are equivalent if they result in the same optimal death benefits, namely, 
D* = D* . However, because the household is maximizing its utility of consumption, it makes more sense 
to say that two risk loadings are equivalent when they result in the same rates of consumption, namely, 
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c\ = Cj for all t > 0. Theorem 5.3 tells us that we achieve that equivalence when the probability of loss 
for the insurance company is the same under the two premium schemes. So, not only is there equivalence 
for the buyers of insurance, the insurers see the two schemes as equivalent from the standpoint of 
probability of loss. 



5.2. Probability of Zero Consumption 

Because we do not restrict the rate of consumption to be non-negative, it is of interest to know 
the probability that consumption reaches zero. We determine this probability for a household buying 
life insurance via a single premium or via continuously paid premium. 

To calculate the probability that consumption reaches zero, we assume that initially the household 
has not bought any life insurance, as in Section 5.1, then purchases the optimal amount of life insurance, 
either D* or D , depending on whether the premium is payable at time zero or continuously until the 
first death, respectively. Correspondingly define tq = mf{t > : c* t < 0}, in which c* is given in (5.3); 
similarly, define r = inf{i > : c* < 0}, in which Z* t is given in (5.6). 

First, assume that the premium is payable at time zero, and calculate the probability P(tq < 
T\ | Cq > 0), in which Cq = r(w — HD*) — — lnk(D*). From Corollary B.3.4 in Musiela and Rutkowski 
(1997), we deduce that the probability density function for r < oo when c* t = Cq + r5t + ^-B t equals 



f T Jt) = T^exp 



■£««+.*>■ 



(5.14) 



with J °° f To (t) dt = e a2rSc o/ m < i because To has a positive probability of taking the value infinity. 
Recall that m=\ f^) 2 , so (i^r) 2 = ^J. It follows that 



POO 

P(r < n I c* > 0) = / f T0 (t) e -^+ x ^ dt = exp 



(5.15) 

The probability P(ri < To < t<i \ Cq > 0) is rather complicated because it depends on whether the 
change in the rate of consumption, Ac* = Ac(D*) from Corollary 3.4, is positive and if it is positive, 
whether it is greater than c*,. Also, note that Ac* depends on whether (x) or (y) dies first, so we write 
Ac* with a — x if (y) dies first and with a — y if (x) dies first. For the sake of brevity, we include only 
an outline of how one calculates P(t! < r < r 2 | c* > 0); then, we write its expression in the Appendix 
for the interested reader. 

If we are given n = t and To > t, then the expression in (5.15), with the drift r8 replaced by m ~ Aa 
and the force of mortality X x + X y replaced by A a , implies that 

P(r <T 2 \T 1 =t, To > t, c*) = e- a < . (5.16) 



From Corollary B.3.4 in Musiela and Rutkowski (1997), we deduce that the joint probability density 
function of c*_ and mino< tt <t c* equals 

. . 2{c-2n + c* Q ) f a 2 r5(n- c%)\ , ( c - 2n + cZ - rdi 
fl( c ' n ) = ,„_^3 7Z ex P 



(^l) 3 i3/2 V rn ) \ ^LVt 
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in which <fi is the probability density function of the standard normal random variable. The domain of 
g is {(c, n) : c > n,n < Cq}. Note that if n = t, then c* = c*_ + Ac*, which might be greater than or 
less than c*_ depending on the sign of Ac*. It follows that 



P(n <r <t 2 \c* >0) = f f f g(c,n)e _a(c+Ac i ) X y e'^^ dcdndt 

J a J Af J c 

JO JAf JC 



(5.17) 



in which we integrate over the two regions in the c-n plane: {(c, n) : c > ( — Ac*) + , c > n, < n < Cq}, 
the first integral with a = x and the second with a = y. Please see equations (A.1)-(A.6) in the 
Appendix for explicit expressions of this probability. 

To compute the probability that consumption reaches zero before the second death add the prob- 
abilities P(to < T\ | Cq > 0) and P(ti < r < t 2 \ > 0). Without loss of generality, suppose that 
Ac* < Ac*. As a function of eg, we have three cases to consider in computing P(r < r 2 | eg > 0), de- 
pending on the relationship of —Ac* < —Ac* with 0. Please see equations (A.7)-(A-9) in the Appendix 
for these three expressions. 

Now, assume that the premium is payable continuously. By analogy with equation (5.15), the 
probability that consumption reaches zero before the first death equals 



P(r < n | Cq > 0) =exp 



(5.18) 



in which Cq = rw — ^ \nk(D*). For the probability that r occurs after the first death and before the 
second, we have the same six cases as detailed in the Appendix. For brevity, we do not write down that 
probability in those cases because the expressions are analogous to those in equations (A.1)-(A.9) with 
Cq, 5, and Ac* replaced by Cq, S, and Ac* = Ac(D*), respectively. 

6. Numerical Example 

We present a numerical example to illustrate some of the results from the previous sections. First, 
set the parameters for the financial and insurance market in this base scenario: 

• The riskless rate of return, or force of interest, is r = 0.02. 

• The drift of the risky asset is ^ = 0.06. 

• The volatility of the risky asset is a = 0.20. 

• The proportional risk loadings are both zero: 6 = 9 = 0. 

Thus, m = i -) 2 = 0.02. Next, set the parameters for the household: 



2 

The force of mortality of (x) is X x = 0.04, so Er x = 25. 

K y = 0.03, so Er y = 33i. 



• The force of mortality of (y) is X y = 0.03, so Et^ = 33,. 



• The income rate of (x) is I x = 2.0, in multiples of $50,000, so the expected lifetime income is 
2(25) = 50 or $2,500,000. 

• The income rate of (y) is I y = 1.5, in multiples of $50,000, so the expected lifetime income is also 
1.5(100/3) = 50 or $2,500,000. 

• The absolute risk aversion of the household is a = 2, which means that the household would pay 
$122.65 to insure against a loss of $10,000 (0.2 of $50,000) with probability 0.01. The expected 
loss is $100, so willingness to pay $122.65 to protect against this random loss seems reasonable. 
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It follows that, per dollar of life insurance, the single premium is H — |, and the continuous premium 
rate is h = 0.07. From (3.5), we compute the optimal death benefit in the case of single premium life 
insurance: D* = 52.38 or $2,619,000. Similarly, from (4.4), we compute the optimal death benefit when 
premium is payable continuously: D = 11.64 or $582,000. Despite the large difference in death benefit 
amounts, we know that these benefits are "equivalent" in the sense of providing the household with the 
same rates of consumption; recall Corollary 5.4. The corresponding probability of loss per policy for 
the insurance company is q = 1 — (7/9) 3 5 = 0.585. 

The change in the rate of consumption at the first death equals Ac* = 0.5476 or $27,380 if (y) dies 
first and Ac* = —0.2024 or —$10,120 if (x) dies first. It is not a surprise that the rate of consumption 
of the household decreases when (x) dies first because of the loss of (x)'s greater income, although the 
death benefit D* = $2,619,000 ameliorates the decrease. The changes in the rate of consumption are 
identical for the case of premium payable continuously, as shown in Corollary 5.4. 

By using the expressions developed in Section 5.2, we conclude from several numerical experiments 
that the probability that consumption reaches zero is so small that one can effectively ignore the anomaly 
of consumption becoming negative. In other words, the cost of the mathematical work involved in 
imposing the constraint that consumption rates be non-negative is not worth the benefit of determining 
the corresponding optimal strategies because they will differ very little from those that we compute in 
this paper. 



7. Summary and Conclusions 

We determined the optimal amount of life insurance for a household of two wage earners. We 
considered the simple case of exponential utility, thereby removing wealth as a factor in buying life 
insurance, while retaining the relationship among life insurance, income, and the probability of dying 
and thereby losing that income. For insurance purchased via a single premium or premium payable 
continuously, we explicitly determined the optimal death benefit; see equations (3.5) and (4.4), respec- 
tively. 

For the reader who wishes to have a handy reference without having to search through the paper 
repeatedly, we reproduce those equations here, along with a list of what the notation means. The 
optimal amount of death benefit, if insurance is charge via a single premium H per dollar of insurance, 
equals 



D* 



max 



1 

ar 



In A x exp al x + 



X x 



+ \ y exp al y + 



A, 



In 



rH 
1 - H 



H 



H 







Here, a is the (constant) absolute risk aversion of the household, X x and X y are the forces of mortality 
for (x) and (y), I x and I y are the rates of income of (x) and (y), and r is the force of interest. The 
optimal amount of death benefit, if premium is payable continuously until the first death at a rate of 
h per dollar of insurance, equals 



D = max 



a(h + r) 



In ! A.,, exp I al x + -^-j + X y exp (al y + -y 



r 







See the Excel spreadsheet at 

http://www.math.lsa.umich.edu/~erhan/OptimalLifeInsurance.xlsx 
for a tool that one can use to calculate either D* or D* . 
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We determined how the optimal death benefit varies with the underlying parameters, and we 
showed that the death benefits are bounded above by annuities based on the maximum income lost; 
see Propositions 3.10 and 4.9. From those propositions, we concluded that the household does not 
over-insure, a desirable outcome of any model. 

In Section 5, we showed that if the premium is determined to target a specific probability of loss 
per policy, then the rates of consumption are identical. Thus, not only is equivalence of consumption 
achieved for the households under the two premium schemes, it is also obtained for the insurance 
company in the sense of equivalence of loss probabilities. 

We welcome the interested reader to modify our model for other common configurations of house- 
holds, as described in Remark 2.1. In future work, we plan to explore optimal life insurance purchasing 
for someone who wishes to bequeath a certain amount of money to a beneficiary. 
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Appendix 

Without loss of generality, suppose that Ac* < Ac*, or equivalently, al y + -y < al x + We 
have six cases to consider, and we provide P(cq) = P(ti < r < r 2 | c*, > 0) for each case. In fact, we 
only have five cases to consider because one can show that Ac* < automatically. However, we include 
case A below for completeness. Technically, case A includes the possible case of Ac* = 0, but one could 
compute the corresponding probability by letting Ac* — > in equation (A.2). 

To simplify the expressions, define 



S=\l{r5Y + ^{\ x + \ y ). 



cr 



A. -Ac* < -Ac; < < c*. 



(2^( S ~ r5 ) + l ) (2 ar5 ~ m ) V / 



B. -Ac: < < -Ac; < c*. 



[0) (^{S-r5) + l){\ar8-m) e \ J 

S V / ^arS — m \ J 



£(S-r5) + a 



(A.2) 
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C. < -Ac* < -Ac; < c*. 



Pc(c* ) = ^-g 

£(S-r*) + a 

+ &(S-r6) + a 
D. -Ac* < < c* < -Ac;. 



Q 2 S A * 



2«r6 — m 



; (4^-«)(-o+ac:) _ 



5 V / ^arS — m \ 



(A3) 



Pd(4) = 



V~" A< e-gV+rSK _ A 

£(S-r8) + i)(\ar6-m) \ J 



(AA) 



E. < -Ac* < c* < -Ac* 

Pe(c* ) = 



A 3 



5 V / iarJ — m 



F. < eg < -Ac* < -Ac;. 

S(£i(S-rd)+a)l y 



(A5) 



(A6) 



Next, we compute p(c*,) = P(ro < T2 | Cq > 0) for the three cases corresponding to the possible 
relationships of - Ac* < — Ac* with 0. 



I. -AcJ < -Ac; < 0. 



p(c* ) = e-^ s -^ c »+P A (c* ), 



in which Pa is given in equation (A.l). 



II. -Ac* < < -Ac;. 



P(4) = 



e -^(s-rS)c* + p D ( c *) ; if o < c* < -Ac;, 
e-^(s-r«)c; + Pb(cS)j if c * > - Ac* 



in which Pd and P# are given in equations (A. 4) and (A. 2), respectively. 



III. < -Ac^ < -a c ;. 



e -%s(S-rS)<i + P F (c*), if < c* < -AcJ, 
p(cg) = <j e -^(s-r«)cS + p B ( c5 ), if _ Ac ; < c* < -Ac;, 
^ e -*£(S-r«)cS + p c ( c *) ) if c * > -a c ;. 



(A7) 



(A8) 



(A9) 
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in which Pp, Pp, and P c are given in equations (A. 6), (A. 5), and (A. 3), respectively. 
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